Height distributions in competitive one-dimensional Kardar-Parisi-Zhang systems 



o 

(N 



in 
o 



■4— » 

o3 



O 

o 



> 
o 
in 
o 

o 

m 



x 

S3 



Tiago J. OliveireQ 

Departamento de Fisica, Universidade Federal de Vigosa, 

(Dated: March 6, 2013) 



36570-000, Vigosa, MG, Brazil 



We study the competitive RSOS-BD model focusing on the validity of the Kardar-Parisi-Zhang 
(KPZ) ansatz h(t) = u^t + (Tt)^x an d the universality of the height distributions (HDs) near the 
point where the model has Edwards-Wilkinson (EW) scaling. Using numerical simulations for long 
times, we show that the system is asymptotically KPZ, as expected, for values of the probability 
of the RSOS component very close to p = p c ~ 0.83. Namely, the growth exponents converge to 
Pkpz = 1/3 and the HDs converge to the GOE Tracy- Widom distribution, however, the convergence 
seems to be faster in the last ones. While the EW-KPZ crossover appears in the roughness scaling 
in a broad range of probabilities p around p c , a Gaussian-GOE crossover is not observed in the 
HDs into the same interval, possibly being restricted to values of p very close to p c . These results 
improve recent ones reported by de Assis et al. [Phys. Rev. E 86, 051607 (2012)], where, based on 
smaller simulations, the KPZ scaling was argued to breakdown in broad range of p around p c . 
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Surface growth is a subject of large interest in Statisti- 
cal Physics, mainly due to pattern formation and emer- 
gency of universality in far from equilibrium systems, 
with applications ranging from biological growth to the 
production of modern thin film devices [J ■ I n stochas- 
tic modeling of the coarse-grained evolution of surfaces, 
the Kardar-Parisi-Zhang (KPZ) equation is one of the 
most important models, which reads 0] 



dh(x, t) 
dt 



vV 2 h+±(Vh) 2 + V (x,t), (1) 



where v accounts for the surface tension, A accounts 
for the growth velocity dependence on the local slope 
and T] is a white noise (with mean null and covariancc 
{rj(x,t)rj(x',t')) = DS(x - x')S(t - t')). Beyond grow- 
ing interfaces, the KPZ class accomplishes several oth- 
ers non-equilibrium stochastic systems, such as directed 
polymers and exclusion processes 0, |3, Q . 

Almost one decade ago, the exact solution of some KPZ 
models in d = 1 + 1 showed that the height evolves asymp- 
totically in time as [f| 0] 



h{t)=v 00 t+{Ttf X , 



(2) 



where Vqq and T arc model dependent parameters, but 
the growth exponent (3 = 1/3 and the asymptotic height 
distributions (HDs) of the random variable x are univer- 
sal. The HDs are associated with the Tracy- Widom dis- 
tributions @ from the random matrix theory of Gaussian 
unitary ensemble (GUE) for curved growth and Gaus- 
sian orthogonal ensemble (GOE) for initial flat interfaces 
S 0] ■ Very recently, the universality of GUE and GOE 
distributions was proven experimentally 0, [13] , numeri- 
cally [nHIl and through the exact solution of the one- 
dimensional KPZ equation with some curved [THEH and 
fiat [HI geometries. 
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Such advances pose interesting new questions in the 
study of growing interfaces. For example, what happens 
to the HDs in competitive KPZ models ? Some works 
have analyzed the skewness of the HDs [13, EH for short 
times, but the first work trying to answer this in the 
light of the Eq. [2] was done very recently by de Assis 
et al. [19| . which studied a competitive growth model 
(RSOS-BD) where particles are deposited with probabil- 
ities p and 1 — p following the restricted solid-on-solid 
(RSOS) and the ballistic deposition (BD) rules, respec- 
tively. Based on scaling exponents from the roughness 
dynamics and in the HDs properties, it was concluded in 
Ref . that there exists three intervals of p with differ- 
ent scaling behaviors: i) Ib ~ [0, 0.75), where the system 
is in KPZ class, with BD-likc surfaces; ii) Ir ~ (0.90, 1], 
where again the KPZ class was found, but with surfaces 
closer to RSOS ones; and Hi) It ~ (0.75,0.90), where 
the system is argued to not be in KPZ class. The growth 
exponent /3 was said to be different from the KPZ one 
and the KPZ ansatz (Eq. ^ not to hold, i.e., it was not 
possible to measure the asymptotic growth velocity 
and the amplitude T. Furthermore, the skewness of the 
HDs was said to perform a continuous transition from a 
positive GOE value in Ib to a negative GOE one in Ir, 
within the It interval. Here, we perform simulations of 
the RSOS-BD model for times one order of magnitude 
larger than the ones studied in Ref. [l{|. Our results 
show that, in It interval, there is a long crossover from 
Edwards- Wilkinson (EW - Eq. [1] with A = 0) to 
KPZ scaling. Thus, indicating that, asymptotically, a 
sudden KPZ-EW-KPZ transition must be expected at 
p = p c m 0.83, where the system is in EW class. An- 
other important result shown here is that the EW-KPZ 
crossover in the roughness scaling seems to be not accom- 
panied by a crossover in the HDs from Gaussian (EW) 
distributions at short times to GOE (KPZ) ones asymp- 
totically. 

It is well known, since a long time, that the competi- 
tion of growth dynamics could produce crossovers in the 



2 



roughness scaling. Generally, in models where particles 
are deposited accordingly to some KPZ rule, with prob- 
ability p, and EW one with probability 1 — p, for small 
values of p there is a crossover in the roughness evolution: 
In early times, the system seems to be in EW class, but 
in a characteristic time t c it crosses over to the KPZ scal- 
ing [2l[ (see Fig. [TJ,). This occurs because the coefficient 
A in the KPZ equation (Eq. Q]) is a continuous func- 
tion of p, which vanishes for p = [22l . [23| . Therefore, 
for small p, A is also small and the laplacian-dominated 
growth produces the apparent initial EW regime. The 
same crossover appears when there is a competition be- 
tween two models in the KPZ class with different signs of 
A in the corresponding KPZ equation. In this case, A = 
for some value p > (see, for example, Ref. [24} )• This 
is the case of the RSOS-BD model: Since A should be a 
continuous function of p interpolating Xbd (positive) and 
Aflsos (negative), there is a value of p (p c « 0.83) where 
A = and the model has the scaling properties of the EW 
class. For any p ^ p c , A ^ and the system must be KPZ 
asymptotically. There are not intermediate universality 
classes continuously interpolating the EW and KPZ ones. 
However, as A — >• (p — > p c ) the crossover time diverges 
with a very large exponent (t c ~ A -4 [HI); thus, it is 
very difficult to observe the expected asymptotic KPZ 
behavior in numerical and/or actual experiments. 

In order to analyze the existence, or possible break- 
down, of the KPZ scaling in the RSOS-BD model close 
to p c , we performed simulations of this model for 0.77 < 
p < 0.89 (the It interval defined above). Out of this 
range of p's, the KPZ scaling was clearly shown in Ref. 
[l9j . Surfaces were grown to times up to t = 2 x 10 5 , 
on a substrate of size L = 2 20 (one time unity corre- 
sponds to L deposition attempts). This large system size 
is important to avoid saturation and finite-size effects in 
our analysis. For each p, averages over 100 samples were 
done. It is important to notice that the RSOS particles 
arc deposited at position i if the constraint (hi—hi±i) < 1 
is satisfied, as in the original definition of the RSOS-BD 
model by da Silva and Moreira [I?}- Using the more re- 
strictive RSOS condition \h% — fo»±i| < 1 the deposition 
of RSOS particles becomes very rare and the EW point 
p c becomes larger than 0.83. 

The surface roughness, or interface width w, is defined 
as the square root of the variance of the HDs. Fig. QJt 
shows the roughness evolution in time for p = 0.81. Af- 
ter an initial transient, typical of BD models, an approx- 
imately EW scaling is found (w ~ t 1 / 4 ). However, for 
longer times, the system crosses over to the expected 
KPZ behavior (w ~ i 1 / 3 ). Obviously, calculating growth 
exponents from linear regressions of a log-log plot, as in 
Fig. into the (nonlinear) crossover region, we could 
wrongly conclude that they arc different from the KPZ 
or EW ones. Indeed, a more reliable measure is provided 
by calculating effective growth exponents (local slopes in 
log-log plots of w x t) as functions of time, as shown in 
Fig. [TJd. For p = 0.83(w p c ) the exponent fluctuates 
around Pew = 1/4 for long times, expected for systems 
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FIG. 1: (color online) (a) Roughness w versus time t for the 
probability p = 0.81. (b) Effective growth exponents /3 e // 
as function of time for several values of p. The dashed lines 
indicates the slopes (8 = 1/4 and /3 = 1/3. 



in EW class. For others probabilities p, the EW-KPZ 
crossover appears as a minimum in the effective expo- 
nents, which approaches 1/4 when p approaches p c . After 
this initial scaling, the exponents converge to the KPZ 
one (Pkpz = 1/3) for all values of p ^ p c studied. As 
discussed before, the crossover time t c , where the KPZ 
scaling begins to dominate, diverges when p —> p c . In 
fact, in Fig. [TJd we notice that the closer to p c , the larger 
is the EW region and the larger are the times needed 
to attain the KPZ regime. For p = 0.81 and p = 0.85 
the exponents did not reach the KPZ even for the large 
times studied, but they are very close to that value and 
still increasing, probably to that. 

From Eq. [2j the non-universal parameter T can be 

obtained from T = [(h 2 ) c / (( X 2 ) c t 2/3 )] ^ where (Y 2 ) c 
is the variance of Y. T becomes constant asymptotically, 
since (h 2 ) c - t 2 / 3 USEE!- This scaling appears for 
long times in the interval It, thus it was concluded in Ref. 
[l9l | that it is not possible to obtain T there. However, 
Fig. UK shows clear plateaus of V for t > 100000 and 
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FIG. 2: (color online) (a) Non-universal parameter estimates: 
r in main plot and Voo in the inset, (b) Skewness of height 
distributions for several values of p in the It interval. The 
dashed lines indicates the expected KPZ GOE values ±0.2935. 



p < 0.81 or p > 0.85. Inside the interval 0.81 < p < 0.85, 
the scaling (h 2 ) c ~ t 2 ^ 3 is not observed in the times 
studied here, as shown in Fig. [lb, possible due to the 
large crossover time. Nevertheless, these results suggest 
that it is always possible to determine T, except at p = p c , 
where the KPZ ansatz (Eq. ^ docs not hold (in other 
words, T = XA 2 jl is identically zero at p c ). 

In the same way, the asymptotic growth velocity (v^ 
in Eq. [5]) could be obtained inside It , but for p very close 
p c only at long times a linear behavior is found in a plot 
of d (h) jdt against t~ 2 l 3 . The inset of Fig. [3^ shows this 
quantity for p = 0.81 and t > 10000. A very good linear 
behavior is found, allowing us to estimate Wqo . We notice 
that, taking the complete time evolution of d (h) /dt, i.e., 
including short times in the graph, a non-linear behavior 
is found, similar to the one shown in Fig. 3a of Rcf. 
fl9j . However, this initial non-linear behavior is due to 
crossover effects and should be discarded in the analysis. 

Finally, in the height distributions (HDs) we still find- 
ing a strong evidence of KPZ behavior in the It region, 



out of p c . Fig. [2Jd shows the skewness evolution for 
several values of p and, at long times, they present a 
good agreement with the Tracy- Widom GOE distribu- 
tion (Sgoe = 0.2935), with exception of values of p very 
close to p c , where S do not converged yet. For p < p c 
it converges to +Sgoe because the system is BD domi- 
nated and this model has a positive A, so that S > 0. On 
the other hand, for p > p c the surface has RSOS-like fea- 
tures and, thus, A < yielding S < 0. For p = 0.83(« p c ) 
the skewness is approximately null, as expected for EW 
systems, whose surfaces have up-down symmetry. We no- 
tice that plotting non- asymptotic values of S (for small 
times) as a function of p, as done in Ref . [l|| , we could 
conclude that there exist a smooth transition between 
+Sqoe and —Sgoe- However, Fig. ^jp shows that in- 
creasing the times where S is taken more abrupt tran- 
sitions will appear, indicating that for t — > oo we must 
find a sudden discontinuous transition at p = p c . 

It is very important to notice in Fig. [2}d that, besides 
the skewness of the HDs to have a minimum for p < p c , 
it seems to not be related to an initial EW regime, in 
contrast to the minimum in effective exponents in Fig. 
[TJ). For p > p c , such minimal (in the absolute value of S) 
does not appear for the studied p's. Therefore, the HDs 
seems do not have a clear crossover from EW to KPZ, i. 
e., they do not cross over from a Gaussian (EW) at small 
times to GOE (KPZ) at large ones. Thus, this suggest 
that the up-down symmetry is broken for any p ^ p c even 
at short times, where the roughness scaling presents an 
approximately EW behavior. However, it could be that 
for p very close p c the crossover appears in HDs, but it is 
very difficult to observe because huge simulations times 
are needed to attain the KPZ regime. Anyway, this is 
different from the roughness scaling, where the crossover 
is present for a broad range of p. 

In conclusion, we have performed large scale simula- 
tions of the competitive RSOS-BD model near the EW 
point p = p c . For several values of p close to p c , 
our results arc consistent with an asymptotic KPZ scal- 
ing, as expected. The growth exponents converge to 
13 = 1/3, but have a minimum associated with an ini- 
tial EW regime. The non-universal parameters and 
r from Eq. [2] could be obtained near p c , but as close to 
p c larger is the time needed to do this. Asymptotically, 
the HDs have a good agreement with the GOE Tracy- 
Widom distribution, but in short times it has absolute 
values of skewness S smaller than Sgoe- All these re- 
sults suggest that there is not a breakdown of the KPZ 
class inside the It interval, in contrast to the proposal 
of Ref. [Hjl, but just a long EW-KPZ crossover that 
makes the KPZ properties to appear at very long times. 
Thus, for t — > oo we must expect a sudden KPZ-EW- 
KPZ transition at p = p c . Another important result 
is that the EW-KPZ crossover in the roughness scaling 
seems do not be accompanied by a counterpart in the 
HDs, where an initial Gaussian (EW, S = 0) regime was 
not observed for the p's studied, suggesting that the up- 
down symmetry is severely broken for all p ^ p c at short 
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(and long) times. This could lead the HDs to converge 
faster than the growth exponents to the asymptotic ones, 
as suggested by Figs. Wp andEb- Therefore, the analysis 
of HDs could be more reliable than scaling exponents to 
determine universality in competitive systems. 
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